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A COMPARISON OF SOLAR CELL DAMAGE BY ALPHA PARTICLES AND PROTONS 


by Albert B. Smith and James W. Blue 
Lewis Research Center 

SUMMARY 

The short circuit current of silicon solar cells illuminated with 2800° K light has 
been measured after the cells received successive increments of proton or alpha particle 
bombardment. The incident particle energies were 10. 5 and 42 MeV respectively, and 
therefore have identical velocity distributions and ranges in the solar cell. The ratio of 
the proton flux to the alpha particle flux to produce the same degradation in short circuit 
current is 3. 8. For particles of equal velocity, the ratio of the Coulomb cross section 
for alpha particles to that of protons to produce silicon recoils of equal energy is 4. If 
the experimental uncertainties are considered, the experimental result is in agreement 
with this ratio, and this agreement implies that the relatively few high energy recoils 
from close alpha particle collisions do not contribute significantly to the cell degradation. 


INTRODUCTION 

This investigation examines the degradation of solar cells from alpha particle bom- 
bardment in order to determine whether the damage rate for these particles is anomalous. 
The modivation for this study arose from the fact that approximately 10 percent of the 
particles of a solar flare at any energy are alpha particles (ref. 1). Consequently, the 
degradation of the solar cells of an interplanetary spacecraft by alpha particles should 
be evaluated. 

While protons and electrons have been used extensively in the degradation of solar 
cells, alpha particles have not. The only other study (ref. 2) was a cursory examination 
with 33-MeV alpha particles. In this work one p on n cell and one rather low efficiency 
(6. 5 percent) n on p cell were irradiated. The conclusions were drawn that the damage 
rates for the two types are similar and that the rate is roughly a factor of 5 more than 
for 19 -MeV protons. 



EXPERIMENTAL ARRANGEMENT 


Since the objective was to examine cells for anomalous damage rates, the simplest 
approach was to make a comparison between alpha and proton damaged cells. The short 
circuit current and open circuit voltage of a cell with approximately 100 milliwatts per 
square centimeter of illumination at 2800° K, were measured after increments of irradi- 
ation. The test conditions were made similar to those of other studies in order to facili- 
tate comparison and prediction of alpha particle degradation for cells not included in this 
study. 

The Lewis Research Center 60-inch cyclotron was the source of particles. The ve- 
locity of cyclotron particles is simply the product of 2ir times the oscillator frequency 
and the exit radius. Since the cyclotron is operated at a fixed frequency and fixed exit 
radius, the alpha particles and protons have identical velocities. This makes the com- 
parison of protons and alpha particles more significant since, when these particles have 
the same velocity, they have the same range. The velocity and energy of 42-MeV alpha 
particles, as they penetrate a piece of silicon, are shown in figure 1. These curves 
were derived from tabulated values (ref. 3) of dE/dx. (Symbols are defined in the ap- 
pendix. ) Radiation defects are initiated by collisions between the primary particle and 
silicon nuclei; the struck nucleus becomes a primary recoil which causes additional dis- 
placements. Most of the collisions between the primary particle and silicon nuclei are 
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(b) Overall view. 

Figure 2. - Apparatus for irradiation of solar cells. 
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not close enough to involve nuclear forces; thus, the Coulomb cross section gives a fair 
approximation to the scattering probability, and curves are shown (fig. 1) for the case of 
the alpha particle to transfer 0. 5- and 2. 0-MeV to silicon nuclei. There is, of course, 
a continuum of recoil energies up to the maximum energy which can be transferred in a 
head-on collision. This value 


^max ^ 


MM 


cr Si w 
\ 2 E " 
m « + M Si ; 


(1) 


is also plotted (fig. 1) as a function of the penetration in the cell. 

A schematic drawing and a photograph of the experimental arrangement are shown in 
figure 2. A stationary cell is used to standardize illumination and is located below the 
line of the beam; the cell to be tested can be positioned so as to intercept the beam, or 
it can be retracted in order to permit the beam to enter a Faraday cup. Figure 3 is a 
photograph of the cell holder and drive unit. The beam passes through a very thin gold 
foil, and a small fraction is scattered into a solid-state detector. The constant of pro- 
portionality between the number of particles hitting the solid-state detector and the num- 
ber collected in the Faraday cup is determined with the cell withdrawn. With the cell to 
be irradiated in the beam, the number of particles hitting the cell is determined by the 
solid-state detector count, the constant of proportionality, and the fraction of the beam 
area that the cell intercepts. Special attention was taken to record the first increments 
of damage at low flux levels. 

The area of the beam was determined by irradiating a piece of X-ray film at the cell 
position; the uniformity was determined by irradiating a copper foil, cutting it into pieces, 
and determining the relative activity of the pieces. The beam could be made uniform to 
better than 5 percent by using quadrupole lenses to defocus the beam. 
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Figure 3. - Motorized cell assembly. 


4 



II 


III 


IV 


<u O' 

1 ^ 

o 

ra — 
"EL o 

- g 

•o i- 
o 

»- <S 

2 


°*Kh 


Region 

It 

v{Q) - 1- 

-f Ed "\r^ 2E d 


Energy of recoil nucleus, Q 


-f- 




_o a 


Figure 4. - Assumed variation of number of displacements with silicon recoil energy. 


THEORY 

The relative number of defects produced by protons and alpha particles will be de- 
termined under the following conditions: 

(1) Minority carriers diffusing from distances greater than 2 Lq (Lq is the initial 
diffusion length) have such a small probability of reaching the junction that a sharp cut- 
off of 2 Lq has been assumed in the calculation. 

(2) The primary particle from the Coulomb interaction produces recoil silicon atoms 
with energies Q, from the displacement energy E^, to Q max - 

(3) The number of displacements per primary recoil v(Q) is given by the Seitz- 
Harrison equation (6) and is shown in figure 4. 

(4) The number of radiation defects affecting cell performance is proportional to the 
number of displaced atoms. 

A primary particle (mass Mj, charge Z incident upon n g ^ silicon atoms per 
cubic centimeter has an energy E(x) after having penetrated some distance x and has 
some differential probability da(E(x),Q) of causing a recoil of energy Q to Q + dQ. 

The number of recoils in this energy range, per square centimeter of surface, for an 
incident particle flux $ is 



x(Q) or 2 Lq 


ct[E(x), Q]dx 


( 2 ) 


The upper limit of the integral is 2Lq when Q < 4|MjM2y/(M^ + MgJ^Efc = 2Lq) and 
otherwise, x(Q), which can be obtained for alpha particles from the curve of Q max ( x ) 

(fig. 1, p. 2). The Rutherford scattering law for an alpha particle on a bare silicon 
nucleus can be substituted for the differential cross section to give 
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Number of particles per initial particle flux, N IQ (cm" 2 )/(MeV" 


( 3 ) 


N(Q)dQ = n gi $ 


2?rZ^e 4 Z 2 


M, 
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x(Q) or 2 Lq 


dx 
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This integral has been evaluated numerically, and the resulting primary recoil spectra 
are shown in figure 5. 

The proton and alpha particle primary recoil spectra are identical in shape and differ 
by just a factor of 4 in number, for Q less than about 1 MeV. Alpha particle bombard- 
ment does differ from protons by producing a continuum of recoil energies to 17 MeV. 

Radiation damage is considered to consist of the silicon atoms displaced by the pri- 
mary recoils. The radiation defects that affect the solar cell behavior appear to be 
(refs. 4 and 5) the result of the radiation produced vacancy or interstitials combining with 
impurity atoms. The number of such defects is assumed to be proportional to the number 
of displaced atoms per square centimeter, which is given by 



Figure 5. - Spectrum of recoil energies for 42-MeV alpha and 10.5 MeV 
proton bombardment of silicon. 


/•^max 

D = J N(QMQ)dQ (4) 

^d 

The dependence of u( Q) is given in 
figure 4 and can conveniently be con- 
sidered in four regions of energy for 
the primary recoil. For region I 
(E^ < Q < 2E^), the primary recoil can 
produce no displacements; consequently, 
v - 1 for the displacement produced by 
the primary particle. For region II 
(2E ( j < Q < E.), the Seitz-Harrison 
equation (ref. 6) is assumed to apply. 

For recoil particles with energy 
greater than the particle is as- 
sumed to lose energy by electron 
ionization and to produce no dis- 
placements until the energy falls to 
Ep at which point it produces i'(E^) 
defects. This energy can be estimated 
in several ways; the one that is taken, 

E. = (M 1 /8m )E (ref. 7), is based up- 
on ionization occurring across the 
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energy gap E^. Region HI extends to the maximum proton recoil energy and region IV 
from this energy to the maximum alpha recoil energy. The number of displacements for 
each of these regions can be obtained in explicit form when it is recognized that the cell 
region 2 Lq is thin enough so that the velocity of the particle in the integral for N(Q) can 
be replaced by its average value. The number of defects per incident alpha particle in 
the four regions becomes 
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( 6 ) 
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( 8 ) 


The number of defects per square centimeter produced in the layer 2 Lq by a 42-MeV 
alpha particle, for the four regions in consecutive order, is 


D tot = ( 8 - 2 + 128 + 9. 4 + 0. 085)$ (9) 

The numer of defects per 10. 5-MeV proton is exactly one-fourth of the alpha particle 
value for regions I, n, and III because there is one-fourth the number of primary recoils 
at any energy below about 1 MeV. For protons, there are no recoils with energies in re- 
gion IV. The number of defects produced by alpha particles in this region is relatively 
small, so it follows that 
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Ratio of short circuit currents, (I /I )' 



Initial particle flux, 3>{cm" 2 )(sec -1 ) 


Figure 6. - Irradiation of several n on p cells. 



Figure 7. * Change in open circuit voltage with proton and alpha bombardment. Resistivity 
10 ohm-centimeters; cells supplied by International Rectifier. 
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This result, when combined with an assumed proportionality between solar cell de- 
fects and the number of displacements, predicts that an alpha particle flux <&q will cause 
the same reduction in cell output as a proton flux of 4<3 >q. 

The low estimate of the damage from region IV for the high-energy primary recoils 
is made under the assumption that has been made for v(Q). It should be recognized that 
there is the possibility of spike phenomena in which the high-energy recoils react with 
many silicon atoms simultaneously. Such effects are seen in fission fragment damage, 
an example of which is seen in the work of Fleischer, Price, and Walker (ref. 8) and, if 
important in this instance, these effects would make the ratio of proton to alpha flux for 
equal damage greater than 4. 


RESULTS AND DISCUSSION 

The results from the testing of a number of n on p cells of materials with different 

resistivities are shown in figure 6. The data are presented as log-log plots of 
2 

CL n/D ■■ 1 against <I> in order to give approximately straight curves of the same slope. 

S y U S 

The alpha particle irradiated cell data scatter in a wider interval than the proton data; 
however, the inclusion of a 1 ohm-centimeter Heliotek cell has no counterpart in the pro- 
ton group. Therefore, the width of the alpha group may not be as wide as shown. A 
comparison of a cell in one group with its counter part in the other group gives an average 
of 3. 8 for the ratio of proton flux to cause the same degradation as an alpha particle flux. 

A test of this equivalence between proton and alpha particle fluxes is shown in fig- 
ure 7, where the degradation in the open circuit voltage for proton and alpha particle bom- 
bardments is shown. In addition, the proton curve has been replotted for flux points at 
<3>p/3. 8. This procedure makes the reconstructed proton curve coincident with the alpha 
particle curve. 

The scatter between cells of the same type makes it questionable whether the differ- 
ence between the expected and observed result, that is, 4 and 3. 8 for the ratio of equiva- 
lent fluxes, is real. Since the experimental ratio is near the predicted value, the high- 
energy primary recoils do not appear to be contributing to cell degradation anymore than 
would be predicted from the Seitz-Harrison theory. In other words, "fission" spike phe- 
nomena from high-energy recoils, if they take place, are not important to cell degrada- 
tion. The lower ratio could be understood in terms of an enhanced annihilation between 
interstitials and vacancies because of the higher density of such defects along the path of 
a primary alpha particle. 

It would be expected that the observed ratio between proton and alpha particle fluxes 
would hold for other energies as long as the comparison is made between protons of 
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energy E and alpha particles of energy 4E, that is, between particles of the same ve- 
locity, for the range of energies where Coulomb scattering is the dominant process. At 
very low energies, for example where the charge exchange is important, the ratio of 
fluxes may be different. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 27, 1966. 
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APPENDIX - SYMBOLS 


D 

number of defects 

°tot 

total number of defects pro- 
duced by protons 

D tot 

total number of defects pro- 
duced by alpha particles 

D“ 

number of defects produced by 
alpha particles in region I 

D a 

u h 

number of defects produced by 
alpha particles in region II 

D a 

number of defects produced by 
alpha particles in region in 

d“ 

U IV 

number of defects produced by 
alpha particles in region IV 

dE/dx 

rate of energy loss, 
MeV/ (g)(cm 2 ) 

E(x) 

energy after penetration of 
x -distance 

E d 

displacement energy, eV 

Eg 

energy gap 

E i 

minimum energy of recoil 
nucleus for which ionization 
is dominant energy loss 
mechanism 

E « 

energy of alpha particle 

e 

charge of electron 

*s 

short circuit current, mA 

I s, 0 

initial short circuit current, 
mA 

L 0 

initial minority carrier diffu- 
sion length, cm 

M Si 

mass of silicon atom 



mass of alpha particle 

Mi 

mass of primary (incoming) 
particle 

m 2 

mass of target (silicon) atom 

m e 

mass of electron 

N 

number of particles 

N(Q)dQ 

2 

number of recoils/cm 

n si 

O 

silicon atoms/cm 

Q 

energy of recoil nucleus, MeV 

^max 

maximum energy of recoil 
nucleus, MeV 

^max 

maximum proton recoil energy 

Q" 

Tnax 

maximum alpha particle recoil 
energy 

R 

range of particle 

V oc 

open circuit voltage, mV 

V oc, 0 

initial open circuit voltage, mV 

V 

velocity, cm/sec 

v(x) 

velocity along x, cm/ sec 

X 

penetration distance 

x(Q) 

penetration distance as a 
function of (Q) 

Z 1 

atomic number of primary 
particle 

Z 2 

atomic number of target nucleus 

V 

number of displacements 

KEj) 

number of displacements per 
primary particle when 
Q = E. 


ll 



u(Q) 

da 

da/ dQ 


number of displacements per 
primary recoil of energy Q 

differential cross section 

differential cross section for 
energy Q, b/MeV 


2 1 

$ particle flux, (cm )(sec ) 

2 1 

$ proton flux, (cm - )(sec ) 

r 

$ a alpha particle flux, 

(cm - ^)(sec - l) 
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